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ABSTRACT 

Aims. We present a model of the global magnetic field in the barred galaxy NGC 1365 based jointly on the large-scale velocity field 
of interstellar gas fitted to H i and CO observations of this galaxy and on mean-field dynamo theory. The aim of the paper is to present 
a detailed quantitative comparison of a galactic dynamo model with independent radio observations. 

Methods. We consider several gas dynamical models, based on two rotation curves. We test a range of nonlinear dynamo models that 
include plausible variations of those parameters that are poorly known from observations. Models for the cosmic ray distribution in 
the galaxy are introduced in order to produce synthetic radio polarization maps that allow direct comparison with those observed at 
AA3.5 and 6.2 cm. 

Results. We show that the dynamo model is robust in that the most important magnetic features are controlled by the relatively well 
established properties of the density distribution and gas velocity field. The optimal agreement between the synthetic polarization 
maps and observations is obtained when a uniform cosmic ray distribution is adopted. These maps are sensitive to the number density 
of thermal ionized gas because of Faraday depolarization effects. Our results are compatible with the observed polarized radio intensity 
and Faraday rotation measure if the degree of ionization is between 0.01 and 0.2 (with respect to the total gas density, rather than to 
the diffuse gas alone). We find some indirect evidence for enhanced turbulence in the regions of strong velocity shear (spiral arms and 
large-scale shocks in the bar) and within 1-2 kpc of the galactic centre. We confirm that magnetic stresses can drive an inflow of gas 
into the inner 1 kpc of the galaxy at a rate of a few M yr~' . 

Conclusions. The dynamo models are successful to some extent in modelling the large scale regular magnetic field in this galaxy. 
Our results demonstrate that dynamo models and synthetic polarization maps can provide information about both the gas dynamical 
models and conditions in the interstellar medium. In particular, it seems that large-scale deviations from energy equipartition (or 
pressure balance) between large-scale magnetic fields and cosmic rays are unavoidable. We demonstrate that the dynamical effects of 
magnetic fields cannot be everywhere ignored in galaxy modelling. 
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1. Introduction 

NGC 1365 is one of the best studied barred galaxies. It has 
been observed in a broad range of wavelengths, including Hi 
(Ondrechen & van der Hulst 1989), molecular gas (Curran et al. 
120011 1, Ha (Lindblad 1999J, and the radio range (Sanqvist et al. 
119951 Beck et al. 2005), in addition to numerous optical and in- 
frared observations (see Lindblad 1999 and references therein). 
Detailed gas dynamical modelling by Lindblad et al. ( 1996} pro- 
vided quantitative models for the gravity and gas velocity fields 
in this galaxy that fit the H i and, to some extent, the CO obser- 
vations. 

The aim of this paper is to add to these efforts by the inclu- 
sion of magnetic fields. The gas dynamical model of the galaxy 
then can be tested against independent radio data, which were 
not included in the construction of the model. Of course, this in- 
volves an additional piece of theory and some further assump- 
tions (concerning, e.g., the applicability of dynamo theory to 
galaxies and uncertainties in some dynamo parameters). Some 
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features of the dynamo theory certainly are not understood well 
enough. However, we demonstrate that gross features of the 
model galactic magnetic field - at least in barred galaxies where 
the shear in the large-scale velocity is the dominant induction 
effect - are rather insensitive to the poorly known details of the 
dynamo system (most importantly, the a-coefficient). Therefore, 
we can plausibly constrain the freedom within the dynamo mod- 
els, and so draw conclusions about the interstellar medium in 
barred galaxies. 

We find fair agreement between radio polarization observa- 
tions and the magnetic field obtained as a solution of the mean- 
field dynamo equations, using velocity and density fields ob- 
tained from gas dynamical simulations, although the distribu- 
tion of polarized intensity is reproduced better than that of po- 
larization angles. Our models also support the idea that inter- 
stellar turbulence is enhanced in the vicinity of dust lanes near 
the bar major axis, and that the energy density of cosmic rays 
can depend only weakly on position in the galaxy, thus devi- 
ating significantly from equipartition with interstellar magnetic 
field. As a result, radio polarization observation and modelling 
of magnetic fields are important ingredients of both the theory 
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and observations of barred galaxies. This work resembles quite 
strongly an earlier study of another barred galaxy, NGC 1097 
(Moss et al. 1200 11 1 . but represents a significant improvement in 
that we now use a dynamical model that specifically models 
NGC 1365, rather than the generic dynamical model adopted for 
NGC 1097. Also, the dynamo model we use here is fully three 
dimensional, whereas that of Moss et al. (f2001 ) used the 'no-z' 
approximation to remove explicit dependence on the vertical co- 
ordinate. Broadly comparable studies have also been published 
by Otmianowska-Mazur et al. (2002), Soida et al. (2006} and 
Vollmer et al. (120061 

2. The observed magnetic structure 

NGC 1365 was observed in total and polarized radio continuum 
with the VLA DnC array at A3.5 cm and A6.2 cm. The full de- 
tails and the maps at 15" and 25" angular resolution are given 
in Beck et al. (2005 ). The total radio intensity (a measure of to- 
tal magnetic field strength and thermal emission) follows well 
the optical bar and the spiral arms. According to the observed 
spectral indices, the thermal fraction is about 20% at A 6.2 cm. 

The polarized emission (Fig.[TJ is strongest in the central re- 
gion and inner bar, but decreases rapidly towards the outer bar. 
There is also significant polarized emission between the bar and 
the spiral arms. No concentration in the spiral arms can be de- 
tected. At /16.2 cm, where the sensitivity is highest, the polar- 
ized emission forms a smooth halo around the bar. The degree 
of polarization is low in the bar and spiral arms, indicating that 
the turbulent magnetic field dominates in the regions of high 
gas density and strong star formation, while the regular field is 
strong between the bar and the spiral arms. At /13.5 cm, most 
of the extended polarized emission outside the bar is lost in the 
noise because of the steep synchrotron spectrum. Furthermore, 
the sensitivity of the VLA to extended structures is reduced for 
scales beyond 3 arcminutes at A3.5 cm, which affects the visibil- 
ity of the large-scale polarized emission in NGC 1365, while at 
A6.2 cm the critical limit is 5 arcminutes and so does not affect 
our observations. 

The peak polarized intensity is 368 mJy per beam at A3. 5 cm 
in the massive dust lane northeast of the centre (see Beck et al. 
2005 ). The fractional polarization is 0.8. At the same position 
the /16.2 cm map reveals a local minimum with polarized inten- 
sity of 150 mJy/beam, corresponding to a fractional polarization 
of only 0.2, which is near the expected contribution from instru- 
mental polarization by the bright nuclear region. This indicates 
that strong depolarization occurs at A6.2 cm in the central region, 
by a factor of at least 4. In the bar and spiral arms the depolar- 
ization factor is 2-3 (Beck et al. 120051 ) . 

Polarized emission can emerge from coherent, regular mag- 
netic fields or from anisotropic random magnetic fields; these 
possibilities can be distinguished with the help of Faraday rota- 
tion measures. In NGC 1097, anisotropic fields dominate in the 
bar region (Beck et al. 2005 ). However, due to the weak polar- 
ized intensity in NGC 1365, the observations available cannot 
provide a large-scale map of Faraday rotation, so that the rela- 
tive contributions of coherent and anisotropic random magnetic 
fields remains unclear. 

3. The model 

3.1. Gas dynamical models of NGC 1365 

We reproduced the gas dynamical model of Lindblad, Lindblad 
& Athanassoula ( 1996 ) using their gravitational potential 'BSM' 



kindly provided to us by P. O. Lindblad. This potential (the 'LLA 
model' in the following) includes the gravitational potentials of 
the disc and spiral arms and was derived from the nonaxisym- 
metric part of the deprojected /-band image. Their best fit pa- 
rameters are Ab ar =1.2 and A sp i ra ] = 0.3 for the relative contri- 
butions of the bar and spiral arms. The model rotation curve fits 
the Hi rotation curve for galactocentric distances r > 120" and 
gives reasonable resonance locations inside this radius. Various 
versions of the LLA model used the bar angular velocity of 
Q p = lSkms^kpc -1 (model BSM) and nkms-'kpc 1 (model 
BSM2), with the corotation radius close to 14kpc in both cases. 

The full gravitational potential of the LLA model is obtained 
from two independent observations: (i) the Hi rotation curve, 
used to fix the total radial mass distribution of the galaxy includ- 
ing dark matter, and (ii) the /-band data, tracing the stellar mass 
distribution, which is only used to derive (after deprojection) the 
amplitude of nonaxisymmetric perturbations in the disc plane. 
The latter cannot be used to derive the rotation curve reliably 
because of the presence of dark matter, and the former also can 
be misleading when the gas flow is significantly nonaxisymmet- 
ric. 

Lindblad et al. ( fl996l > adopted 20 Mpc (1" = 97 pc) for the 
distance of NGC 1365, but we adjusted the model to a distance 
of 18.6 Mpc (1" = 90 pc) (Lindblad [T999l 

Isothermal gas dynamical models were calculated using the 
code ZEUS 2D, published by Stone & Norman JT992] ), and w e 
found a close match to the model of Lindblad et al. Q996). 
However we did not attempt to take into account the warp in 
the outer disc, as we are mostly interested in the inner region. 
Our basic models, illustrated in Fig. |2j),c have the bar angu- 
lar velocity Q p = 16.16kms _1 kpc _I and the corotation radius 
at R c = 15.5 kpc; we also considered a model (Fig. [2^) with 
Q p = 17 km s~' kpc 1 and R c = 16.3 kpc. The angular velocity 
of the spiral pattern is taken to be equal to that of the bar. For 
reasons explained below in Section [4] the resulting gas density 
in the bar region was too low to reproduce the observed mag- 
netic field within the dynamo model. The gas density in the LLA 
model can be argued to be underestimated inside the corotation 
radius because the rotation curve used had poor resolution, and 
underestimates the depth of the potential well. We derived our 
basic model from the LLA model, by replacing the rotation curve 
used by Lindblad et al. (1996) with the more recent CO rotation 
curve of Sofue et al. i 19991 ). This modified model was much bet- 
ter able to reproduce the observed magnetic field, while remain- 
ing in agreement with the overall morphology of the molecular 
gas distribution. A significant difference is that there is more ma- 
terial in the central regions when Sofue's rotation curve is used. 
The rotation curve used here is shown in Fig. [3^, with the posi- 
tions of resonances illustrated in Fig. [3}). 

We also studied the dependence of the gas dynamics and 
magnetic field on the sound speed adopted in the isothermal gas 
model. This parameter is uncertain in our models for several rea- 
sons. Englmaier & Gerhard (1997 ) showed that the large-scale 
gas distribution in isothermal gas flow models of barred galax- 
ies can depend on the sound speed, even if the pressure forces 
are negligible. Since the position, and even existence, of shocks 
depends on Mach number, the global gas flow configuration can 
change as a result of a relatively small change in the speed of 
sound. Different parts of the multi-phase interstellar medium 
(ISM) may not follow the same global gas flow. Different nu- 
merical methods have been shown to represent different aspects 
of the ISM with varying success. Sticky particle methods, for 
example, model better the clumpy ISM, while grid-based meth- 
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Fig. 1. The polarized intensity contours and magnetic vectors of the polarized radio emission at the wavelengths /13.5 cm (left 
hand panel) and /16.2 cm (right hand panel) (both smoothed to a resolution 25"; the beam size is shown in the lower right of 
each panel), superimposed onto an ESO optical image of NGC 1365, kindly provided by P. O. Lindblad. The contour levels are 
1, 2, 3,4, 6, 8, 12, .. . times 30/iJy/beam at A3. 5 cm and 40 ^Jy /beam at A6.2 cm; the r.m.s. noise is 15 /iJy/beam at A3. 5 cm and 
14yuJy/beam at /16.2 cm. 
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Fig. 2. The model gas density with superimposed velocity vectors in the reference frame corotating with the bar, in gas dynamical 
models based on (a) the rotation curve of the LLA model with c s = lOkms 1 (left hand panel), and the rotation curve of Sofue et 
al. ((T999) with (b) c s — 10 km s~' (middle panel) and (c) c s = 30 km s _1 (right hand panel), with c s the sound speed. Shades of grey 
represent the logarithm of gas density (darker shades corresponding to larger values), with each shade corresponding to the same 
density in each panel. Note the smaller density contrast in the bar region in the model with higher speed of sound (panel c). 



ods give a better description of the shocks and the smooth gas 
component. 

The global magnetic field depends on the gas flow via Eqs ([TJ 
and <j2j; however, it is not a priori clear which component of 
the ISM carries the magnetic field and, therefore, what is the 
appropriate sound speed of the gas. We have considered models 
with the speed of sound equal to 10 and 30 km s~' (see Sect. 15.41 ). 



Our magnetic field model also relies on the gas density ob- 
tained from gas dynamical simulations together with the velocity 
field; this is discussed in Sect. 13.21 - see Eq. (01). 
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Table 1. Parameters of models discussed in the text, as defined in Sect. 13.21 In all the models, the angular speed of the bar is 
O p = 16.16kms _I kpc 1 with the corotation radius at 15.5 kpc. 
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3.2. The dynamo model 

Dynamo models, specifically simple mean-field turbulent dy- 
namos, are remarkably successful in explaining the observed 
features of galactic magnetic fields (see, e.g., Ruzmaikin et al. 
[19881 Beck et al. IT9961 Widrow l2TX)2l for reviews). Despite 
the fact that the nonlinear behaviour of turbulent dynamos is 
still controversial, mean-field models provide a remarkably reli- 
able empirical description of large-scale (regular) galactic mag- 
netic fields (Shukurov 2004). Fortunately, dynamo solutions for 
galaxies are quite insensitive to those parameters that are poorly 
known, such as the form of the a-effect and even, to a lesser 
extent, the turbulent magnetic diffusivity. This is especially true 
of models for barred galaxies where large-scale velocity shear 
plays a dominant role in determining magnetic field structure 
(Moss et al. [1998a 2001); then the primary role of the ar-effect 
is to maintain the field against decay. 

Our model can be regarded as a development of the dynamo 
model of Moss et al. ( 120011 1, used to model the large-scale mag- 
netic field in a generic barred galaxy. We now introduce fur- 
ther elaborations required to reproduce the basic features of the 
global magnetic pattern in NGC 1365. We solve the standard 
mean field dynamo equation for the large-scale, regular mag- 
netic field B 



dB i 

dt 1 



aB 



\VrjxB 



rjWxB), 



(1) 



in three spatial dimensions, using Cartesian coordinates (x, y, z), 
where x and y are horizontal dimensions, and the disc midplane 
is at z — 0. Here a and rj are the turbulent transport coefficients 
responsible for the a-effect and turbulent magnetic diffusion, re- 
spectively, u is the large-scale velocity field, and the term with 
V77 allows for the turbulent diamagnetism associated with the 
spatial variation of the turbulent diffusivity (Roberts & Soward 
119751 1. In our standard case, our computational domain covers 
the region -L < (x,y) < L, < z < ah — z max , where a is the 
domain's aspect ratio. We take a mesh of size n x x n y x «,, with 
uniform spacing in the horizontal directions and also, separately, 
vertically. The maximum resolution readily available to us was 
n x = n y = 200, n z = 31, and in order to resolve satisfactorily the 
solutions we took L = 15 kpc and a = 0.12, so z max = 1.8 kpc. 
(Thus we study only the inner part of this unusually large barred 
galaxy.) The total thickness of gas layer that hosts the large-scale 
magnetic field is taken as 2h = 0.9 kpc, compatible with the 
thickness of the diffuse warm gas in the Milky Way. Our proce- 
dure is to time-step the x and y components of Eq. (Q}, and then 
to use the condition V ■ B = to update B z . We restrict ourselves 
to solutions of even (quadrupolar) parity with respect to the disc 
plane z = 0, and so the latter step is straightforward, given that 
B, — at z — 0. This is the same procedure used in the three- 



dimensional galactic dynamo models described in Moss i 19971 1, 
except that cylindrical polar coordinates were used there. 

In Eq. (Q]l, a parameterizes the dynamo action of the inter- 
stellar turbulence, and 77 is the turbulent magnetic diffusivity. We 
assume both of these quantities to be scalars (rather than tensors) 
and, in order to obtain a steady state with saturated dynamo ac- 
tion, introduce a simple a-quenching nonlinearity into the prob- 
lem, writing 



O'O 



1 + %B l /B, 



fi 2 = \np{r)v\ 



eq 

O(r) 

a o = a f -^—f(z) , 



(2) 



(3) 



with 



sm(nz/h) . 



|z| < h/2 . 



f(z) = 



[cosh(2|z|//z- l) 2 ] 1 sgnz , \z\ > h/2 



Here Qo is a typical value of Q, Z? eq is the magnetic field strength 
corresponding to equipartition between magnetic and turbulent 
kinetic energies, and a* is a constant, which we can adjust. Quite 
arbitrarily, we adopt Qo = Q. at r = 3 kpc, and Eq. ([3]) shows that 
a* is the maximum value of a at this radius. Thus we are as- 
suming that the large-scale magnetic field significantly reduces 
the (T-effect when its energy density approaches that of the tur- 
bulence; the constant £ is introduced to suggest formally some 
of the uncertainty about the details of this feedback. The de- 
pendence of a on height, defined by f(z), is implicitly odd with 
respect to the midplane, with |or| increasing with \z\ from at 
z = to a maximum at \z\ = h/2, and then decreasing to zero as 
\z\ — > 00 (remembering that we only explicitly model the region 
z > 0). Because of the symmetry of Eqs ([TJi and (f2]), if B is a 
solution, then -B is also a solution. 

We take ^ = 0(1), assuming that there is no catastrophic 
ff-quenching (Brandenburg & Subramanian 2005). The models 
were computed with f = 1, and the field strength then scales as 
£-'/ 2 . The turbulent speed that enters B eq is taken to be equal 
to the speed of sound as adopted in the gas dynamical model. 
The gas density p(x, y, 0) is taken from the gas dynamical model 
described in Sect. 13.11 We extend this away from z — by writing 



p(x,y,z) = 



p(x,y, 0) 
cosh(|z|//z) 



(4) 



The magnitude of the gas density is relatively unimportant in our 
model (where the Lorentz force is not included into the Navier- 
Stokes equation) as it affects only the magnitude of the magnetic 
field in the steady state, via Eq. (0, but not its spatial distribu- 
tion. The only aspect where the magnitude of gas density plays a 
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Fig. 3. (a): The rotation curves used in the paper: that from 
Lindblad et al. (1 19961 1 (solid; as in Fig. [2^), and one more con- 
sistent with more recent CO observations (Sofue et al. 1999) 
(dashed; as in Fig. |2]3,c). The plot assumes the distance of 
NGC 1365 to be 20Mpc as in Lindblad et al. ( fl99"6l ). The ra- 
dius of corotation is R c * 14 kpc. (b): The linear resonance 
diagram for the rotation curves shown in (a) with the same 
line style. From bottom to top: £2 - k/2, Q, and O + k/2 in 
units of km s -1 kpc 1 . The resonances are located at the in- 
tersections with the horizontal lines corresponding to Q p = 
16.16 km s -1 kpc~' (solid) and 17kms _1 kpc _1 (dashed). The 
small scale structure in the Q. + k/2 curves is an artefact of plot- 
ting. 

role is the Faraday depolarization and, hence, the modelled dis- 
tribution of polarized intensity. This effect is, however, relatively 
weak at A = 3-6 cm and it is plausible that other depolarization 
effects (e.g., Faraday dispersion) are more important in the real 
galaxy. The gas density in our model is shown in Fig. [4] 

The gas velocity in the plane z = 0, u(x,y, 0), is also taken 
from the gas dynamical model. For convenience, we split this 
into rotational and non-circular parts, 

u(x, y, 0) = £l(r)r<f> + v(x, y, 0), (5) 

respectively, where r — (x 2 + y 2 ) 1 ^ 2 is axial distance. 

We then introduced two significant modifications. We found 
that, in the gas dynamical model, Q(r) increases very rapidly 
towards the rotation axis (very approximately, as 1/r). The gas 
dynamical model model appears to handle this feature satisfacto- 
rily, but it causes significant numerical problems for the dynamo 



code at attainable numerical resolution. Thus Q was softened 
by introducing an explicit parabolic profile within a radius of 
2. 1 kpc, with the maximum of Q truncated to 1 10 km s _1 kpc -1 
(as compared to 1730 km s~' kpc -1 at r = 0.013 kpc, the small- 
est distance from the axis in the gas dynamical model used). 
This modification can be expected to reduce the magnetic field 
strength in regions close to the galactic centre, but as this region 
is not well resolved by the radio observations, we cannot in any 
case make a comparison between these and the computed mag- 
netic field. 

Further, we continued the velocity field above the disc by in- 
troducing z-dependence into the horizontal velocity components 



u(x,y,z) = 



u(x,y,Q) 



cosh(|z|/1.2kpc) 



(6) 



and u z = everywhere. 

In order to model a galaxy surrounded by near-vacuum, we 
allow the magnetic diffusivity to become large high in the halo 
(Sokoloff & Shukurov H990l . 



'?() 



1 



1 + 071-D 



1 - exp - 



\z\-h 
1.5 kpc 



< h 



, \z\>h 



where 770 and 771 are constants; thus 77 = 770 near the disc midplane 
and 77 — > 770/71 in the halo region (|z| > h). We adopted a nominal 
rji = 2 - larger values led to numerical difficulties. A conven- 
tional value of 770 is 10 26 cm 2 s _1 ; however, we also considered 
models with values larger than that - see Table Q] In order to 
reproduce polarized radio maps of NGC 1365 in sufficient de- 
tail, we had to introduce further spatial variation in 77. Following 
Moss et al. (2001), we have assumed that the turbulent diffu- 
sivity is enhanced by the shear of the nonaxisymmetric velocity 
according to 



S 



Sr 



du. 



dy 



8Uy 

dx 



where S max is the maximum value of S . The effect of f n + 
is, firstly, to broaden magnetic structures near the spiral arms, 
and, secondly, to reduce the central peak of magnetic field. The 
values of /„ adopted are shown in Table Q] We did not consider 
a similar enhancement in a as Moss et al. (2001) found it to be 
unimportant. The values of f l; that were sufficient to produce re- 
alistic magnetic fields in spiral arm were still too small to reduce 
the central maximum of magnetic field to an acceptable level. 
Therefore, we introduced an additional enhancement of 77 in the 
central part of the galaxy, multiplying 770 by q n exp(-r 2 /2r 2 ); the 
values of q n and r n are given in Table[T]for each model studied. 

Clearly, we have made a number of rather arbitrary choices, 
in particular when extending the two dimensional gas dynam- 
ical model into three dimensions. Our overall impression from 
a substantial number of numerical experiments is that the over- 
all nature of our results does not depend very strongly on these 
choices. 

Atz = Zmux, and on x,y = ±L, the boundary conditions are 
B x = B y = 0. On z = 0, dBJdz = dB y /dz = 0, B, = 0, and so 
the integration of V • B =0 gives the values of B, on the other 
boundaries. These are conservative boundary conditions on B x 
and By, in that they will increase the field gradients and thus 
raise the threshold for dynamo action to occur. 

We nondimensionalize the problem in terms of the length 
L — 15 kpc, time /i 2 /77o and magnetic field B eq . Given that the 
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velocity field, including the angular velocity, is given by the dy- 
namical model, the only free dynamo parameter is a*; the corre- 
sponding dimensionless parameter is 

aji 

Ra = —, (7) 

m 

where a, is defined in Eq. (0). The dynamo action prevents mag- 
netic field from decay for values of R a exceeding about 1 for 
rjo = 10 26 cm 2 s _I ; the critical value of R a increases roughly pro- 
portionally tO J](). 

Henceforth, we will use dimensionless variables, unless ex- 
plicitly otherwise stated; the units of gas number density and 
magnetic field strength are 44 cnT 3 and 30 yuG, respectively. 

4. Results 

The gas dynamic and dynamo models described above together 
yield the gas density and the distribution of the large-scale mag- 
netic field in the galaxy. The distribution of magnetic field in 
the galaxy plane resulting from Model 2 (introduced in TableQ]), 
which we argue below to be our best model, is shown in Fig. [4] 
Given the distribution of the cosmic rays, we can now construct 
synthetic radio observables in order to assess the quality of the 
model. We have computed synthetic radio polarization maps at 
wavelengths of 3.5 cm and 6.2 cm using the dynamo generated 
magnetic field and the gas density, and compared them with 
the observed radio maps. Details of this procedure are given in 
Appendix [Al Since we do not model turbulent magnetic fields, 
we are unable to calculate the total radio intensity and to estimate 
the degree of polarization from the model. 

We considered several models for the number density of cos- 
mic rays, n cr , which we discuss in Sect. 15. II For all but one of the 
dynamo models listed in Table[T]we find that the simplest possi- 
ble choice, n cr = const, provides the best fit to the observed data, 
regardless of the other quantities adopted. The larger value of r n 
in Model 1 produces a relatively weak magnetic field through- 
out a large central region compared with that at the ends of the 
bar. In order to fit the observed central peaks of polarized in- 
tensity, P, an implausible non-uniform distribution of n cr would 
be required in this model. Specifically, the cosmic ray distribu- 
tion required to reconcile this model with observations would 
have a high peak within 3 kpc of the centre where magnetic field 
strength is minimum. For the other dynamo models, any plausi- 
ble non-uniform distribution of « cl produces too strong a central 
maximum of P relative to all other structures. In particular, the 
polarized intensity in the spiral arms is almost lost in models 
with non-uniform n cr , being far weaker than that within 1-2 kpc 
of the centre. Since we have truncated the angular velocity at 
r < 2.1 kpc, the untruncated differential rotation would lead to 
an even stronger discrepancy. 

Our synthetic maps do not include any depolarization effects 
due to random magnetic fields (see Burn 1966; Sokoloff et al. 
I1998K although they allow fully for depolarization by the regu- 
lar magnetic fields (differential Faraday rotation and beam depo- 
larization). In order to include Faraday depolarization effects due 
to the large-scale magnetic field, we assumed a nominal constant 
ionization fraction of X = n e /n = 0.1, corresponding to a ther- 
mal electron density of 0. 1 of the total gas density obtained from 
the gas dynamical simulations as described in Sect. 13.21 Guided 
by analogy with the Milky Way, where the average total gas den- 
sity is 1 cirT 3 whereas the thermal electron density is 0.03 cm -3 , 
a smaller value of X might be appropriate. We show results for X 
close to this value in Sect. 14. II In Sect. 15.51 we discuss the effect 
of variations in X and argue that 0.01 5 X 5 0.2. 
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x (kpc) 

Fig. 4. Energy density contours and vectors of the regular mag- 
netic field B from Model 2 (see Table Q}, both at z = 0, are 
shown together with gas number density represented with shades 
of grey. The contours shown correspond to approximately 0.1, 
0.6 and 3.0 times the r.m.s. value; the length of the vectors is 
proportional to B 1 . The scale bar at the top of the frame refers to 
the gas number density in the units of hydrogen atoms per cm 3 . 



We used two main techniques to compare the synthetic maps 
with observations and therefore to select the optimal magnetic 
field model. We chose to use the Ab.2 cm map of polarized in- 
tensity in the analysis since it has the best signal-to-noise ra- 
tio. All model data, including synthetic radio maps, have been 
smoothed (in terms of the Stokes parameters Q and U) to match 
the resolution of the observations. In Sect. 14.21 we compare the 
distributions of polarized intensity on cuts along various paths 
in the plane of the sky. In Sect. 14.31 we analyse the difference 
between the computed and observed polarized intensities in two 
dimensions. In addition, we compare the orientations of the mag- 
netic B-vectors obtained from the observed and synthetic Stokes 
parameters (Sect. [43). 

To rotate the model galaxy to the position of NGC 1365 in 
the sky, we took the inclination angle i = 46° and the position 
angle of the galaxy's major axis (i.e the intersection of the sky 
plane and the galaxy plane) PA = 222°, which are those assumed 
in obtaining the rotation curve for our (favoured) gas dynamical 
model. Results are quite sensitive to these values, and it is pos- 
sible that a reappraisal could result in noticeable changes. 

4.1. Synthetic polarization maps 

Overall, Model 2 (specified in Table [TJ appears to provide the 
best fit to the observed polarization map; Model 4 is only slightly 
worse - see Sect. I4.2I Contours of B 2 shown in Fig. [4] indicate 
that the regular magnetic field is stronger in the bar region where 
gas density is large, and outside the regions of high density in the 
spiral arms. There are magnetic features apparently unrelated to 
the density distribution [e.g., those passing through the positions 
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Fig. 5. A synthetic map of polarized synchrotron intensity (con- 
tours) and polarization planes at A6.2 cm, resulting from Model 2 
(see Table[TJ assuming that n CI = const, are shown superimposed 
on the optical image of the galaxy NGC 1365 (shown in only 
a few shades of grey for clarity). The synthetic map has been 
smoothed to the resolution of 25" to match that of the observed 
map shown in Fig. Q] The contour levels shown are approxi- 
mately (1,3,6, 12,32) x Pmax/45, where P max is the maximum 
of P in the synthetic map. Dashed lines show the position of cuts 
discussed in Sect. 14.21 



(x,y) * (-5, 8), (5, -8)]; they are presumably formed by a lo- 
cally enhanced velocity shear. The magnetic field has a deep 
minimum within the bar, mainly produced by the density defi- 
ciency in that region. Other important features clearly visible in 
Fig. |4] are the magnetic field enhancements in the dust lane re- 
gion, where magnetic field is amplified by both compression and 
shear, and the prominent central peak. 

The synthetic polarization map for this model is shown in 
Fig. [5] This can be compared directly with the observed map in 
the right-hand-panel of Fig. [TJ the maps (and all other maps we 
show) are at a similar scale to facilitate the comparison; we make 
this comparison more quantitatively in Sect. 14.31 Our models 
have a high degree of symmetry, whereas the 'real' NGC 1365 
is only approximately symmetric; since the observed map looks 
more regular on the eastern side, we shall mostly refer to that 
part of the galaxy unless stated otherwise. Despite the differ- 
ence in symmetry, there is broad agreement between these two 
maps; for example, both have a deep minimum of P near the 
bar's major axis where gas density is low, and both have the 
magnetic spiral arms displaced from the gaseous ones (although 
both magnetic arms are displaced to larger radii in the synthetic 
map, only one arm is so displaced in the observed map). The 
minimum of the synthetic P in the bar (corresponding also to a 
minimum of magnetic field within the bar, as seen in Fig. [4), 
is broader than of the observations (see Sect. 14. 2\ . The rea- 
son for this is the very low gas density in this region, lead- 
ing to weaker magnetic fields via Eq. (|2). This feature is fur- 
ther discussed in Sect. [6] where we argue that the gas dynam- 




-100 100 

distance from centre (arcsec) 




distance from centre (arcsec) 



Fig. 6. Cuts, at position angle -31° passing through the galac- 
tic centre (left to right in the plots corresponds moving from 
south-east to north-west in the sky), through polarized intensity 
maps at A6.2 cm smoothed to HPBW = 25", for (a) the observed 
map, and synthetic maps from (b) Model 2 and (c) Model 4, both 
for n a - = const. In panels (b) and (c), the synthetic profiles for 
/16.2 cm and /13.5 cm are shown solid and dotted, respectively; 
the difference is due to Faraday and beam depolarization for the 
assumed ionization degree X = 0.1. The units of P are as in 
Fig. [Tjfor (a) and arbitrary in (b) and (c), but adjusted to fit a 
similar range. The dotted profiles for A3. 5 cm with X — 0.1 also 
correspond to P at A6.2 cm with X = 0.032. 



ical model underestimates significantly the amount of molecu- 
lar gas in the bar region. Synthetic P is large both to the north 
and south of the bar major axis. In particular, the model repro- 
duces a maximum of P upstream of the bar major axis, centred 
in the /16.2 cm map of Fig.[TJat (RA = 03h 33 min 40 sec, Dec = 
-36°09'00"). These maxima apparently arise from slightly en- 
hanced velocity shear (that locally amplifies magnetic field) 
rather than from local density maxima. We also note maxima 
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Fig. 7. As in Fig. [6] but at position angle 0° (left to right is south Fig. 8. As in Fig. [6] but at position angle -90° (left to right is 



to north in the sky). 



east to west in the sky). 



of P near the ends of the bar and the beginning of the spiral 
arms, at (RA = 03h 33 min45 sec, Dec = -36°08'15") and 
(RA = 03h 33 min 28 sec, Dec = -36°08'30"). We note t hat the 
observed total emission (not shown here; see Beck et al. 2005) 
is related to gas density in a rather straightforward manner being 
correlated with the gas density. The fact that this is not the case 
with the polarized intensity (as seen in both observed and syn- 
thetic maps) confirms that the observed regular magnetic field is 
not frozen into the gas, apparently being affected by the dynamo 
action. 

4.2. Cuts through polarization maps 

We found that comparisons can be usefully quantified and de- 
tailed using the cuts in the sky plane mentioned above. We show 
cuts only through the map at /16.2 cm because this map has 
higher signal-to-noise ratio and includes the large-scale emis- 
sion fully. However, depolarization is significant at this wave- 



length (see Sect. |2]i and has to be taken into account when com- 
paring the model and observations. We use cuts through the 
centre of the galaxy at position angles PA = 0°, -90° and 
-31°, where PA is measured counterclockwise from the north 
as shown in Fig. [5] (The angle -31° is chosen so that the cut 
goes through the spiral arms; this corresponds roughly to a di- 
agonal in the computational frame of Fig. |4]) The positions of 
these cuts are shown in Fig. [5] The synthetic P has been nor- 
malized to make the mean difference between that and the ob- 
served P approximately zero. We have superimposed another 
profile from cuts through synthetic maps which represents both 
P at /13.5 cm with our favoured value of X = 0.1, and also 
at /16.2 cm with X = 0.032 (incidentally, this is close to the 
mean ionization degree of the warm diffuse gas in the Milky 
Way). The coincidence of these two cuts is due to equivalent 
Faraday depolarization, which depends directly on the quantity 
iff(z) °c A X J nB\\dz along the line of sight towards an observer 
(see Appendix [A}. Here B\\ is the line of sight field component 
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and n e = Xn, therefore A 2 X = const identifies equivalence in 
depolarization. We see that this value is about the same in both 
cases (i.e. 6.2 2 x 0.032 * 3.5 2 x 0.1). The difference between the 
polarization for these two possibilities is then just a /l-dependent 
scale factor. 

The cuts are presented in Figs |6j [7] and [8] for the best-fit 
Model 2 and also for Model 4. The latter model has the back- 
ground turbulent magnetic diffusivity rjo enhanced by a factor 
of 2. This leads to a significantly smoother, less structured dis- 
tribution of P. Thus, comparison of Models 2 and 4 allows us 
to suggest that the effective turbulent magnetic diffusivity in 
the interstellar gas of barred galaxies is, on average, close to 
rjo = 10 26 cm 2 s~'. This value is typical of spiral galaxies in 
general and is that obtained if the turbulent speed v t is close to 
10 km s" 1 and the turbulent scale is about / = 0. 1 kpc; 770 - \lv t . 

Our model neglects depolarization due to random magnetic 
fields which can reduce the value of P in the central parts more 
strongly than in the outer galaxy and therefore affect the relative 
height of the central peaks in Figs [6H5] Depolarization due to 
internal Faraday dispersion reduces the degree of polarization to 

1 -e- s 

P=Po ~S~' (8) 

where S = 2cr 2 M /i 4 with cr 2 M = 2C\{b 2 ){n 2 )dL the variance of 
the Faraday rotation measure. Here C\ is the dimensional con- 
stant appearing in the definition of the Faraday rotation mea- 
sure (see Appendix [A}, b is the turbulent magnetic field, angu- 
lar brackets denote averaging (the fluctuations in magnetic field 
and thermal electron density are assumed to be uncorrected), 
d is the turbulent scale and L is the path length (Sokoloff et al. 
1998). The best available estimate of the random magnetic field 
in the central region of NGC 1365, b 40 //G, follows from 
the total synchrotron intensity assuming equipartition between 
cosmic rays and magnetic fields (see however Sect. |5.1| for a dis- 
cussion of the validity of this assumption). For n e = 0.03 cm" 3 , 
d — 0.1 kpc, L — 1 kpc and A = 6.2 cm, we then obtain S - 6, 
implying that this mechanism can depolarize the central peak 
significantly, giving p/po - 0.2. Since the height of the sec- 
ondary peak should also be affected by depolarization, albeit to 
a lesser extent, we expect that the ratio of the two peaks will 
be reduced by a factor smaller than five. We note, however, that 
this estimate is uncertain since the number density of thermal 
electrons, their filling factor, turbulent scale and other parame- 
ters are not known well enough. An alternative is to assess the 
importance of this depolarization effect by comparing polarized 
intensities at A6.2 cm and A3.5 cm. The ratio of the central peak 
to the secondary ones at A3.5 cm is about 6-8, as opposed to 2-3 
at A6.2 cm. The difference can be attributed to Faraday depolar- 
ization (by both regular and random magnetic fields). Assuming 
that Faraday depolarization at A3. 5 cm is negligible, we conclude 
that it can reduce the degree of polarization at A6.2 cm by a fac- 
tor as large as 4, which is consistent with the analytical estimate. 
We conclude that the relative height of the central peak in the 
synthetic cuts of Figs |6H5] would be reduced by Faraday disper- 
sion, although this is difficult to estimate accurately. 

Given the above uncertainties in the amount of depolariza- 
tion, all three cuts for Model 2 are similar to those observed. In 
particular, the relative heights of the peaks in P and, more im- 
portantly, the positions of both maxima and minima are remark- 
ably realistic. The characteristic feature of this model is that rj 
is further enhanced by a factor of q n - 3 in the inner region 
of NGC 1365, r 5 3 kpc. This enhancement can be due to a 
higher rate of star formation, and hence more hot gas, with a 



correspondingly higher speed of sound, which would allow the 
turbulent speed to be larger than elsewhere. 

The Model 2 cut at PA = -31° (Fig. [6), which passes through 
the spiral arms, shows an encouraging agreement with observa- 
tions. For example, B has a maximum slightly outside the north- 
ern arm in both this model and the real galaxy. However, the out- 
ermost maxima produced by the spiral arms are slightly too far 
away from the centre in the model. As illustrated in Fig.[6jb), the 
relative heights of the peaks at A6.2 cm are significantly affected 
by Faraday rotation even for X = 0.1, where they clearly differ 
by more than just a scale factor between A6.2 cm and A3.5 cm. 

The cut at PA = (Pig. O exhibits similar degree of agree- 
ment with the observations. The main deficiency of the model 
here is the too narrow distribution of P (the magnetic structure 
of the model is too poor outside the bar) and the minimum is too 
deep near the centre of the cut. 

The cut at PA = -90° in the synthetic map, shown in Fig. [8] 
has a central maximum that is too narrow (or off-centre minima 
that are too broad). This difference results in the deep minima 
in the difference parameter 6 discussed in Sect. 14.31 The sharp 
minimum in the observed cut near the centre is a result of beam 
depolarization; it occurs in the synthetic cuts as well, but is re- 
moved by smoothing. 

Model 2 seems to be almost optimal. The model could 
be fine tuned by changing 770 and r n within the ranges (1- 
2) x 10 26 cm 2 s" 1 and 1.5-3, respectively. For example, the sec- 
ondary peaks in the PA = cut decrease in strength in Model 4. 
Further, increasing n cl by a factor of 2 within the central 1 .5 kpc 
would make the central peak higher. However, we have not made 
such post hoc adjustments. 

4.3. The difference maps 

To obtain a global comparison of the models and observations, 
we produced maps of the difference between the observed and 
synthetic polarization at /16.2 cm, with the synthetic polariza- 
tion scaled to make the mean difference approximately zero; this 
measure was further normalized by dividing the difference by the 
appropriately normalized noise level of the observed map giving 

( 1 •4-P/-f > niax)model — (^ > /^ > max)obs /r . x 

6 = 7 — 75 — ^ ■ (9) 

Thus, all comparisons were performed pointwise after their re- 
duction to the common resolution 25" - this is quite a stringent 
test of the model. The result is shown in Fig. |9]for Model 2. 

Since the models - unlike the real galaxy - possess perfect 
symmetry, the difference can hardly be uniformly small: a per- 
fect fit in one half of the galaxy would produce significant sys- 
tematic discrepancy in the other half. With this caveat, the dif- 
ference map shows an acceptable global agreement of the model 
with observations, in that it does not show much of the basic 
morphological elements of the galaxy. The normalized relative 
difference is about 6-14 in four spots observed to the east, south 
and north-west of the galactic centre, indicating that synthetic 
polarized intensity is too small upstream of the dust lanes and 
at two positions at the inner edge of the western spiral arm. 
Otherwise, \S\ ^ 4 across the whole field of view. Given the 
limited scope of our model (e.g., it does not include any turbu- 
lent magnetic fields which can produce polarized radio emission 
where they are anisotropic), we consider this degree of agree- 
ment to be acceptable. We discuss in Sect. |5.1| a cosmic ray dis- 
tribution that would provide perfect fit of Model 2 to observa- 
tions. 
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Fig. 9. The difference 5 between normalized synthetic and observed polarization maps at /16.2 cm, as defined in Eq. (0, superim- 
posed on the optical image of NGC 1365. The contour spacing is 2, with the zero contour shown solid, negative values of 6 dashed, 
and positive, dotted. 



4.4. Faraday rotation 

We can use polarized intensity (as in the comparisons above) to 
probe the distribution of the large-scale magnetic field strength, 
and also to deduce the orientation of the magnetic field in the 
plane of the sky (via polarization vectors). However, knowl- 
edge of this quantity does not determine the field direction. The 
Faraday rotation measure RM is sensitive to the direction of the 
magnetic field, but the observed RM map is very patchy because 
of the lower signal-to-noise ratio at A3. 5 cm. Therefore, we used 
RM data only to establish a minimum acceptable degree of gas 
ionization. 

4.5. Magnetic field structure 

An analysis of the observed global magnetic structure in 
NGC 1365 that is sensitive to the direction of magnetic field was 
performed by Beck et al. (2005) by fitting the polarization angles 
obtained from multi-frequency observations. This analysis pro- 
vides the large-scale magnetic field expanded into Fourier series 
in the azimuthal angle. Their results indicate the presence of a 
significant component with the azimuthal wave number m = 1 at 
almost all distances from the galactic centre. However, our un- 
derlying gas dynamical model has even symmetry in azimuth, 



so that modes with odd values of m do not occur in the mod- 
elled magnetic field. The contribution of the m — 1 mode to 
those Fourier expansions is more important than just produc- 
ing the overall asymmetry. In particular, superposition of vari- 
ous azimuthal modes produces local magnetic features at kilo- 
parsec scale which are lost if only even modes are retained in 
the observed structure to facilitate comparison with the model. 
Therefore, we did not find it useful to compare the modelled and 
observed magnetic structures in this manner. (The presence of 
unmodelled odd-ra structure was also a feature of our study of 
NGC 1097 in Moss et al. |200Tl ) 

We instead compare directly the orientation of the magnetic 
field vectors in the observed and synthetic polarization maps. 
Comparison of two-dimensional vector fields is difficult. We 
could approach this by taking cuts through maps of the magnetic 
field orientation angles, as was done with the polarized intensity. 
However, a small shift in a feature such as a shock front can re- 
sult in drastic differences between any such cuts made parallel 
to the front. 

In Fig. [10] we show the orientation of both the synthetic and 
observed magnetic field vectors obtained directly from the cor- 
responding Stokes parameters; points below 3 times the r.m.s. 
noise level are neglected in the observed maps. Agreement be- 
tween model and observations is reasonable in the top left and 
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Fig. 10. Orientations of the B-vectors of polarized emission ob- 
tained from the observed (white) and synthetic (black, Model 
2) Stokes parameters at /16.2 cm. The vector lengths are pro- 
portional to yP, For the observed B-vectors, positions where 
polarized intensity is weaker than three times the noise level, 
P < 3o>, are neglected. 



bottom right quadrants near the bar (and partly further out), 
whereas the differences are quite large in the other regions. The 
difference is especially large near beginning of the spiral arms. 
The mean value of the difference between the observed and mod- 
elled polarization angles is 33°, and its standard deviation is 25°. 
For comparison, the errors in the observed polarization angle 
range from 2 to 10 degrees. 

The overall difference is that the model polarization vectors 
are arranged into a more elliptical pattern around the bar than 
the observed ones, which exhibit a more circular configuration 
(and have almost constant pitch angles). It seems that the non- 
axisymmetric distortion due to the bar is weaker in the observed 
magnetic field than in the model. This could be because the mag- 
netic field is coupled to a warm gas component which has less 
response to the bar's potential than cold gas and stars. We made 
a comparison similar to that in Fig. [K)]but for Model 6, where 
the speed of sound is 30kms _1 (see Table[T|i. The resulting gas 
dynamical model illustrated in Fig. [2] (right hand panel) has a 
more uniform density distribution and weaker deviations from 
axial symmetry. The improvement in the magnetic pattern was 
only marginal, and so the reason for this discrepancy remains 
unclear. 



5. Sensitivity to parameters, and implications of the 
dynamo models 

In this section, we discuss how synthetic radio maps are affected 
by various changes in our model. This allows us to infer useful 
information about the interstellar medium in the galaxy. 
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Fig. 11. Cuts through radio maps for various models of cosmic 
ray distribution as defined in Eq. dTOb ; the position angle of the 
cut is indicated in the upper right corner of each panel: (a) -3 1 °, 
(b) 0, and (c) -90°. The solid lines are for Model 2, i.e., with 
n cr = const, and are identical to those in Figs. |6lb)-[8lb); the 
other two cuts in each panel are for model (ii) « cr =oc yj, 
dashed; and model (iii) n cl = p~ l , dotted. 



5. 1. Distribution of cosmic rays 

In order to calculate the synchrotron intensity, we need to specify 
the number density of cosmic ray electrons. There are no direct 
measurements of this quantity in external galaxies, and there are 
no sufficiently detailed theories that might supply it. Therefore, 
empirical models for the cosmic ray number density n CI were 
studied with 

(i) n a = const, 

(ii) n cr oc V7, (10) 

(iii) n cr oc p~ l 7 
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where / is the observed total radio intensity and p is the percent- 
age polarization of synchrotron emission at short wavelengths; 
models (ii) and (iii) are motivated below. Thus, we make an 
attempt to relate n cr to various observable quantities. The first 
model is the simplest possible, and it attributes all the varia- 
tion in synchrotron intensity to that of the magnetic field. The 
motivation for model (ii) is that, if cosmic rays are in equipar- 
tition with the total magnetic field, w cr oc B 2 oV then the total in- 
tensity of synchrotron emission is roughly proportional to n 2 .. 
Model (iii) relies again on the idea of equipartition between 
cosmic rays and magnetic fields, but now we use expression 
p ~ p B 2 /B 2 ot (Burn [19661 Sokoloff et al. [1998l > to estimate n ct 
as « C r K B 2 0t =! B 2 po/p, where po w 0.75 and B can be taken 
from our dynamo model. This rough estimate neglects any de- 
polarization effects and any contribution of anisotropic random 
magnetic fields to polarized intensity (Sokoloff et al. 119981 . 

As illustrated in Fig.QT| the major and universal effect of any 
plausible non-uniform distribution of « cl is to enhance the cen- 
tral maximum in P, so that the peaks in the outer parts become 
relatively insignificant. Most importantly, any signature of the 
spiral arms almost disappears. If the synthetic polarized inten- 
sity is rescaled to fit that observed in the spiral arms, the central 
peak becomes unacceptably broad and high. Using n cr oc B 2 /p 
instead of p as in model (iii) does not improve the situation. 
Similarly, it does not help if we use n cr oc B 2 ot — B 2 + 4npv 2 , 
assuming that the random magnetic field is in equipartition with 
turbulent energy. Since models (ii) and (iii) involve the assump- 
tion of energy equipartition (or pressure balance) between cos- 
mic rays and magnetic fields, we conclude that our results do not 
suggest this type of relation between cosmic rays and magnetic 
fields at large scales. 

We can think of several plausible explanations for this per- 
haps surprising result. It may be that some fraction of the po- 
larized emission in the outer galaxy is produced by anisotropic 
magnetic fields which are not modelled. We do not consider 
this to be a very plausible option as this would require that the 
anisotropy is larger in the outer galaxy and in the spiral arms, 
rather than in the region of the central peak. Meanwhile, velocity 
shear, which might produce the anisotropy, is stronger in the in- 
ner bar region. More plausibly, cosmic ray diffusion makes their 
distribution smoother than that of the magnetic field. With the 
cosmic ray diffusivity of K =s 10 29 cm 2 s _1 and the confinement 
time t =i 10 6 yr, their distribution would be rather homogeneous 
at scales (Kt) 1 ' 2 - 1 kpc. We note, however, that our model sug- 
gests that cosmic ray distribution is almost uniform at scales of 
order 10 kpc. 

It cannot be excluded that the synthetic polarization maps 
exaggerate the relative height of the central peak because they 
neglect Faraday depolarization due to random magnetic fields, 
namely the internal Faraday dispersion discussed in Sect. 14.21 
As follows from Eq. ®, and discussion following it, it is not 
implausible that this effect can reduce the relative height of the 
central peak by a factor of five or somewhat less. As can be seen 
from Fig. QT| the ratio of the central to the secondary peaks at a 
distance of about 50" from the centre is about 20 or more for the 
non-uniform distributions of cosmic rays, whereas the observed 
ratio is about 2-3. Given the uncertainty of any estimates of the 
amount of depolarization, we cannot exclude that models with 
a non-uniform distribution of cosmic rays could be reconciled 
with observations. 

We show in Fig.Q~2]the ratio of polarized intensity observed 
at A6.2 cm to the integral along the line of sight J B\ ds, with B± 
the component of the modelled magnetic field in the sky plane. 




x (arcsec) 



Fig. 12. A map of the ratio of polarized intensity observed at 
A6.2 cm to the synthetic P obtained from Model 2 with n cr = 
const. This ratio can be used to assess the variation of n CT re- 
quired to achieve a perfect fit to observations (neglecting any 
anisotropy of the turbulent magnetic fields and any contribution 
of a hypothetical galactic halo to the polarized emission). The 
contour at level unity is shown solid, the other contours are at 
levels 3", with negative and positive integer values of «; contours 
above (below) unity are shown dashed (dotted). 



The latter is proportional to the synthetic polarized intensity ob- 
tained for a uniform cosmic ray distribution. If our magnetic field 
model were perfect, the above ratio would show the variation of 
cosmic rays across the galactic image. We note that the value 
of the ratio varies remarkably little in the bar region; the most 
prominent variations arise from the local peaks of the observed 
polarized intensity that are also prominent in Fig. [9] Figure [T2l 
confirms that the variation of cosmic ray energy density within 
the galaxy is rather weak and consists of a large-scale, smooth 
variation with contours of a shape similar to that of gas density 
and other tracers in the bar, and perhaps with a few local max- 
ima. 

5.2. The effects of the turbulent magnetic diffusivity 

Models 2 and 4 yield very similar magnetic field distributions 
(see Fig. [6]), even though the background turbulent magnetic dif- 
fusivity T]o in Model 4 is more than twice as large as in Model 2 
(see Table [TJ. The main effect of enhanced magnetic diffusivity 
is to make the secondary peaks of P in Model 4 less prominent, 
even with n cr = const. Model 4 could be reconciled with ob- 
servations if n CI were enhanced in the outer bar regions and in 
the spiral arms and/or reduced in the central part. Unless this is 
the case, and given that Model 2 agrees with observations bet- 
ter than Model 4, we conclude that our models support a value 
T] Q <, 10 26 cm 2 s- 1 in NGC 1365. 

One of the effects of turbulent magnetic diffusion (and dy- 
namo action) is to produce a misalignment between magnetic 
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Fig. 13. The alignment of the regular magnetic and velocity 
fields (in the non-rotating frame) illustrated with the contours 
of cos^ = \u ■ B\/uB at the midplane for Model 2, projected on 
to the optical image in the plane of the sky. The levels shown are 
0.4 (dotted), 0.7 (dashed) and 0.95 (solid). 



field and velocity in a shearing flow. Since the velocity shear is 
strong everywhere in the bar region and near the spiral arms, we 
expect that magnetic and velocity fields would be tightly aligned 
(in the corotating frame) if magnetic field were frozen into the 
flow (Moss et al. 120011 

The degree of alignment between the model magnetic and 
velocity fields in NGC 1365 is illustrated in Fig. [13] where we 
note that the angle between the two vectors exceeds 20° almost 
everywhere in the bar. The misalignment is also significant near 
the spiral arms. 

The local enhancements of turbulent magnetic diffusivity, by 
a factor of 2-3, in the dust lanes and near the galactic centre in- 
troduced in Sect. |3.2| are important in our model as they allow us 
to avoid excessively large strength of the large-scale magnetic 
field produced by extreme velocity shear in those regions. As 
argued by Moss et al. (2001), such a local enhancement of in- 
terstellar turbulence may be associated with instabilities of the 
shear flow. 



5.3. The role of dynamo action 

Model 3 has the same parameters as Model 2, except that R a = 0. 
Thus, this model includes the same effects of rotation and ve- 
locity shear as Model 2, but without any dynamo action. In 
this model, the magnetic field decays on a timescale of about 
0.6 Gyr, after an initial accelerated decay near the galactic cen- 
tre. However, the magnetic structure in the outer parts of the 
galaxy is remarkably similar to that with R a + 0. We deduce 
that the magnetic structure we have obtained does not depend 
strongly on details of the poorly known a-effect, but rather is 
controlled by the large-scale velocity field, which is known much 
more reliably. The role of the a-effect is just to maintain the mag- 
netic field against decay, which is enhanced by the strong shear 
typical of barred galaxies. This situation is similar to that found 
when modelling another barred galaxy, NGC 1097 (Moss et al. 

USD. 



5.4. The effect of the speed of sound 

Models 2 and 6 have different values for the speed of sound (10 
and 30 km s _1 , respectively). The higher speed of sound results 
in less structure in the velocity and density fields (even though 
we neglected to include an equivalent increase of the turbulent 
magnetic diffusivity by a factor of 3, which would be neces- 
sary for strict consistency). The cuts for Model 6 have relatively 
weaker features in the outer parts of the galaxy. This result ap- 
pears less acceptable, and we deduce that 10km s _1 is a more 
favourable value for the speed of sound of the gas phase, to 
which the regular magnetic field is coupled. 



5.5. Dependence on the gas ionization fraction 

For X ^ 0.2, the synthetic cuts of P show a much greater dis- 
agreement with the observed cuts than when X = 0.1. For exam- 
ple, for the cuts in Fig. [8] the effect of increasing X (and therefore 
increasing depolarization) is to broaden the minima in P about 
+50", and increase the ratios of the central maximum to outer 
maxima (at about ±100"). 

When we calculate synthetic RM for values of X 55 0.01, 
the range of values obtained does not match the observed range 
±600 rad rrT 2 . For example, at X = 0.0 1 , the maximum synthetic 
RM is about 350 rad irT 2 . 

Given that a constant ionization fraction is appropriate, our 
models suggest that 0.01 S5 X % 0.2 to be a plausible range of 
values. 



6. Discussion and conclusions 

We have constructed a three dimensional dynamo model for 
NGC 1365, with the rotation curve, non-circular velocities and 
gas density taken from a dynamical model for this particular 
galaxy. Thus, although we have taken a similar approach as in 
our earlier studies of other specific barred galaxies (IC 4214, 
Moss, Rautiainen & Salo [T99%l NGC 1097, Moss et al. |200TT l, 
for NGC 1097 we adopted a generic dynamical model as in- 
put, whereas here we have been able to use a bespoke model. 
We have tried to make a much more detailed comparison be- 
tween observations and model predictions than previously - but 
see also a recent paper by Vollmer et al. (I2006I I. 

Of course, we have been restricted to using a mean field dy- 
namo model - for this sort of study there is really no plausible 
alternative currently available. Our modelling (and that of the 
earlier papers, cited above) has demonstrated that when mod- 
elling galaxies with strong non-circular velocities the role of the 
mean field a-coefficient is primarily to offset the inevitable dif- 
fusive decay of the field, and thus allow a steady state with fields 
of order equipartition strength to be maintained. The major de- 
terminant of the field structure is the non-circular velocity field 
(Sect. [53] and Moss et al. |1998al |2001| l, and the main features 
can be expected to persist for plausible field maintenance mech- 
anisms. We must further bear in mind the other limitations of 
the modelling, including the restriction to the inner part of the 
galaxy, r < 15kpc, which means that boundary effects may in- 
fluence results near this radius. 

Our main conclusions are as follows. We see no evidence for 
a variation in n cr as strong as in B 1 . This may imply that equipar- 
tition between cosmic rays and the regular magnetic field is not 
maintained even at global scales. The discrepancy between our 
crude models involving the equipartition assumption and obser- 
vations can be significantly reduced if Faraday depolarization 
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due to turbulent magnetic fields is taken into account. We have 
discussed this result further in Sect. 15.11 

The strongest deviations of the synthetic polarized intensity 
from that observed occur in the bar region, just upstream of the 
dust lanes. The reason for the low synthetic polarized intensity 
is the small value of magnetic field strength there in the model. 
More precisely, the depth of the minima in synthetic and ob- 
served P are similar but they occur at somewhat different posi- 
tions, and those in the synthetic map are broader (see Fig. [8]). 
However, the relative heights of the maxima in P are reproduced 
quite successfully. In this sense, the agreement is better than 
might be inferred from Fig. [9] The reason for the difference is 
the deep and broad minimum in the gas density in those regions. 
We believe this to be a shortcoming of the gas dynamical model, 
which was fitted to incomplete CO data. In particular, CO ob- 
servations of NGC 1097 (Crosthwaite 120011 ) do not show the 
minima of the density in the bar region to be as deep as in the 
model of NGC 1365 used here. 

Rather surprisingly, the agreement between the orientations 
of the model magnetic field and the B-vectors derived from the 
observed polarization vectors is not as good as that between the 
model and synthetic PI distributions. We have discussed this in 
Sect. 14.51 The reason for this difference remains unclear. 

Our preferred model relies on the galactic rotation curve 
and gas density distribution different from those suggested by 
Lindblad et al. |1996t the rotation curve used is that resulting 
from CO observations (Sofue at al. |1999l >. Our results are com- 
patible with the observed distributions of polarized synchrotron 
intensity and the magnitude of the Faraday rotation measure for 
the number densities of ionized diffuse gas of order 0. 16 cnT 3 at 
a distance of order 5 kpc from the centre along the bar's minor 
axis and 0.21 crrT 3 in the spiral arms. With the gas dynamical 
model used here, this corresponds to the mean ionization frac- 
tion of 0.01-0.2. 

Our models confirm that magnetic field strengths in the inner 
bar region can be strong enough to drive mass inflow at a rate of 
several solar masses a year (see also Moss, Shukurov & Sokoloff 
120001 Moss et al. |200Tl Beck et al. |2005l . Thus, in these strongly 
barred, strongly magnetic galaxies, it becomes necessary to in- 
clude the dynamical effects of magnetic fields in order to repro- 
duce all features of the gas flow. It follows that self-consistent 
magnetohydrodynamic modelling of barred galaxies is required. 

Keeping in mind our dynamical model is incomplete, at least 
in that azimuthal structure corresponding to odd modes is omit- 
ted, our general conclusion is that mean field dynamo models 
are reasonably successful in modelling magnetic fields in this 
barred galaxy. Moreover, such models can also provide informa- 
tion about both the gas dynamical modelling process and con- 
ditions in the interstellar medium (see also the models for the 
'normal' spiral galaxy M31 in Moss et al. [1998b ). A robust con- 
clusion is that, contrary to widely held opinions, dynamical ef- 
fects of magnetic fields cannot be everywhere ignored in galaxy 
modelling. 



Appendix A: Synthetic radio maps 

Synthetic radio maps have been obtained by computing the 
Stokes parameters Q and U from the magnetic field obtained 
from the dynamo simulations. The magnetic field was rotated to 
the same position in the sky as the galaxy NGC 1365, and Q 
and U at a given wavelength were obtained by integration along 



the line of sight (where the z-direction here points towards an 
observer at infinity): 

Q = C f e(r)cos[2^(r)] dz, 



U = Cj e(r) sin[2(A(r)] dz, 

with allowance for Faraday rotation by the thermal ionized in 
the local polarization angle tfr (see, for example, Sokoloff et al. 
(1998 ) and references therein for more details). Here C is a di- 
mensional constant whose specific value is inessential here, and 
e is the synchrotron emissivity. We assume that the synchrotron 
spectral index is a constant, q - -1, so that e oc n a B\. We use 
models of n cr from Sect. 15.11 

Depolarization due to differential Faraday rotation is in- 
cluded as we take 

n t {r)B\\(r)dz, 

where A is the wavelength of the emission and C\ = 
0.81radirT 2 cm 3 /iG _1 pc~' is a dimensional constant. The num- 
ber density of thermal electrons n e was obtained from the to- 
tal gas number density n as given by the gas dynamical model, 
assuming a constant ionization degree of X. (So we take n e = 
Xp sas /mu with inn the proton mass.) 

The synthetic Stokes parameters are then given by convolv- 
ing the map of raw parameters with a Gaussian beam in sky 
plane. 
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